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The formation of black hole low-mass X-ray binaries: through 

case B or case C mass transfer? 
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Abstract. The formation of low- mass X-ray binaries containing a rather massive (M ^ 7Mq) black hole is 
problematic because in most recent stellar evolutionary calculations the immediate progenitors of these black 
holes (Wolf-Rayet stars) lose so much mass via their stellar wind that their final masses are well below the 
observed black hole masses. We discuss the recently proposed solution that these binaries are formed through case 
C mass transfer (i.e. mass transfer after core helium burning is completed), avoiding a long Wolf-Rayet phase 
and thus significant mass loss. We show that only some of the currently available models for the evolution of 
massive stars allow this formation channel. We also investigate the effect of the downward revised Wolf-Rayet 
mass-loss rate as is suggested by observations, and conclude that in that case Wolf-Rayet stars end their lives 
with significantly higher masses than previously found and may be able to form a black holes. 
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1. Introduction 



evolutionary calculations resulted in the conclusion that 
even massive single stars might end their evolution as rela- 



In low -mass X-ray binaries a neutron star or a black hole ^-y^^y ipw-mass objects when they explode (^challer et al 



-mass (M ;^ IMq) companion, a sec- 
relarively wide bi- 



accretes from a low- 
nario to form such srars begins with 
nary of a massive star and a low-mass companion. When 
the massive star becomes a giant, mass transfer is unsta- 
ble and a common-envelope forms in which the companion 
spirals down towards the core of the giant, leaving a close 
binary consisting of t he helium core of the giant and the 
low- mass companion ( van den Heuvel 1983| ). The helium 
star explodes in a supernova and depending on the (core) 
mass of the helium star, a neutron star or black hole is 



1992|; Meynet et al. 1994; Woosley et al. 1990 and are 



formed. With the discovery of A0620-00 (Eylesetal. 1975; 



Elvis ernl. 19751) Ahd Lhc dcKimiiiiaLloii of the maijtj fuiic 



tion of 3.18 dMcCUiiLock fc RcmlUard 1980[ ), Lhtj exlsLciice 
of the class of black hole low-mass X-ray binaries was es- 
tablished. Currently we know 6 to 8 such systems depend- 



thus unable to produce the observed black holes (see also 
Kalogera 1999 ). For massive stars in close binaries, which 
lose their hydrogen envelopes due to mass transfer early in 
their evolution the situation is even worse; the most recent 
calculations predict masses of helium stars as they explode 
as low as 3 M(i), almost in dependent of their initial mass 
( IWellstein fc Langer 1999| ). 

In this article we first discuss the formation of black 
hole low-mass X-ray binaries through case C evolution as 
suggested by Brown et al. (1999) and Wellstein & Langei 
(19991): mass transfer starting after core-helium burning 



ing on the membership criteria (Charles 1998; Bailyn et al 
199S). An evolutionary scenario for these objects is given 
in |dc Kool et al. (19871 ). 



To make a black hole, the initial mass of the primary 
must exceed a criti cal value, wh ich currently is believed to 
be around 20 Mf7^( Fryer 1999 ). However, large mass-loss 
rates for massive stars and W olf-Rayet stars have been 
inferred from observations (e.g. de Jager et al. 198q ) and 
are found from the comparison of Wolf-Rayet models with 
these observations ( Langer 1989a ). Applying these rates to 
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has been completed (Kippenhahn & Weigert 1967). In 
this case a long-duration Wolf-Rayet phase in which the 
star loses a lot of mass is avoided (Sect. 0). Then we dis- 
cuss the most recently observed mass-loss rates for Wolf- 
Rayet stars and the implication of lower mass- loss rates on 
the final helium-star masses of exploding stars in binaries 
(Sect. H). At the end we discuss uncertainties and possible 
alternatives for the formation of black hole low-mass X-ray 
binaries (Sect. |j) and end with our conclusions (Sect. ||). 

2. Case C mass transfer 

It has been suggested that case C mass transfer could be 
invoked to avoid a long-duration Wolf-Rayet phase in the 
evolution of the massive star, in order that this star does 
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not lose too much mass and still is able to form a massive 



black hole ( [Brown et al. 1999| ; [Wellstein fc Langer 1999| ). 

The occurrence of case C mass transfer depends on 
the radius evolution of massive stars. For supergiants the 
radius of the star is not very well defined, since the outer 
layers of the giant envelope are extremely dilute. However, 
the best we can do is use the calculated values of the 
radii of giants. We also neglect the interaction between 
the wind of the massive star and the companion which 
may influence the separation of the two stars. 

We calculate the initial separation with which a binary 
should start in o rder to undergo case C mass transfer as 
follows (see also Portcgics Zwart et al. 1997 ). The sep- 
aration at the moment the Roche-lobe overflow (RLOF) 
starts is given by 

orlof = — (1) 

where R is the radius of the star and tl is the dimension- 
less Roche-lobe radius (the ratio of the R oche-lobe radius 
and the binary separation). We use the Eggleton (1983 ) 
equation for tl. For mass ratios between 10 and 50, the 
value of tl is between about 0.6 and 0.7. During the evo- 
lution the star loses mass and the separation increases 
according to 

«'-|^. (2) 

where M denotes the total mass of the binary. So to start 
Roche-lobe overflow at time t when the star has a radius 
R{t), the initial separation is given by 



— orlofIO 
R.{t) M{t) 
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possible anymore. The two other panels in Fi g, p] show the 



same, but for the stellar evolution models of Hurley et al 
(2000[]middle) and [Heger et al. (200C| , bottom). For these 
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We now compute the separations at which massive 
stars fill their Roche lobes as function of initial mass and 
initial separation and determine whether the mass trans- 
fer is case B or case C. In Fi g. [l| (top) we show th is for the 
evolutio nary calculations of Schaller et al. (1992 , see also 
Fig. 4 of Kalogera fc Wcbbink 1998 ). For a star of initially 
15 MQcase C mass transfer occurs for initial separation 
between 1000 and 1320 R0. For a 20 M0Star, these lim- 
its are 1300 and 1550 Rr;^. For a 25 Mfi^star case C is not 



Mzams (Msun) 
Fig. 1. Initial separations Oj for which case B and case 

C mass transfer occur as fu nction of ZAMS mas s, for a 1 

M0 companion. Top for the Schaller et al. (1994 ) models, 

mi ddle for the Hurley et al. (200C ) models, bottom for 

the Heger et al. (2000| ) models! 



3. Case B mass transfer 

A different way to avoid too much mass loss may be the 
fact that observed mass loss rates (which are the basis 
for the mass-loss rates used in the evolutionary calcula- 



tions) are revised downward (Hamann & Koesterke 1998 
[Nugis fc Lamers 2000 ) , which may make it possible to pre- 



models case C is not possible for stars more massive than 
around 19 M0. A recent estimate of the number of black 
hole low-mass X-ray binaries that can form through the 
narrow case C interval of the Schaller models shows that 
even such a narrow interval might be enough to explain 



the whole Galactic population (Brown et al. 2001). 

We conclude that since case C evolution depends 
strongly on the radius evolution of massive stars which is 
very uncertain, it seems possible but is not certain whether 
black hole low-mass X-ray binaries can be formed in this 
way. 



vent helium stars in binaries to lose so much mass they 
no longer can become black holes. In a recent paper with 
drastically lower mass-loss rates, derived from one partic- 
ularly well-studied object and extrapolated, final masses 
over 20 M0for the m ost massive helium stars are found 
( Cherepashchuk 2001). As shown by Kalogera (1999] ) the 
helium stars that were the progenitors of the black holes 
in binaries cannot have lost more than half of their initial 
mass. This includes both mass loss in the stellar wind and 
in the supernova explosion. 

A recent compilation of observed mass-loss rates for 
Wolf-Rayet stars is made by Nugis fc Lamers (2000 ). In 
Fig. g we show these inferred mass-loss rates for WN and 
WC/WO stars (excluding the hydrogen rich Wolf-Rayet 
stars). We overplotted mass- loss rates for WN and WC 
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Fig. 2. Mass- loss rates for Wolf-Rayet stars as observed 
(triangles for WN and circles for WC/WO stars, from 
Nugis fc Lamers 2000| , where we excluded the hydrogen 
rich Wolf-Rayet stars) and various relations used for evo- 
lutionary calculations. The solid lines are the relations 
assumed by [Woosley et al. (1995| ) for WC (upper) and 
WN (lower) stars. The dashed lines are the ones used by 
Wellstein fc Langer (1999| ), where we converted the mass 
loss - luminosity relations to mass loss - mass relations us- 
ing the mass - luminosity relation of Langer (1989b). The 



upper dashed line is for their standard case, the lower for 
their reduced mass loss case. The dash-dotted line is a 
rectangular least square fit to all points (see also the text 
and Eq. (|)). 



stars as used by Woosley et al. (1995) as the solid lines 



The mass-loss rates used recently by Wellstein fc Langer 
(1999|]]are shown as the dashed lines, where we used the 
luminosity - mass relation as given by Langer (1989b ) to 
convert the mass loss - luminosity relation used by these 
authors, to a mass loss - mass relation. The top dashed 



Fig. 3. Final helium star masses as function of the initial 
helium star mass with the mass-loss rates according to 
Woosley et al. (1995|, solid line, some of their results are 



potted as solid triangles) and Wellstein & Langer (1999 



dashed line) assuming a helium star lifetime as given by 
Woosley et al. (1995| ) . A selection of their results is plotted 
as the open triangles. The numbers at the top give an 
estimate of the ZAMS mass of the progenitor of the helium 
star. 



the final helium stars mass Mf can be computed from the 
initial mass M\ and the helium star lifetime (r) from 



Mi = [Af/-"-|-(a-l)/cr 



1/(1-") 



(6) 



fine is tlicii' HLaiidai'd case, the boLLom a reduced maay-loHa 
rate. Which Lhey used Lo account foi' the lowei' observed 
mass-loss rates. 

The most recently determined mass-loss rates thus sug- 
gest that the rates used by Wellstein & Langer (1999) are 
still too high. We will investigate the effect of using a 
lower mass-loss rate law, which is shown in the figure as 
the dash-dotted line and is given by 



As a check of our calculations we show in Fig. H the fi- 
nal masses that we obtain using the top two lines (dashed 
and solid ) from Fig. p| and th e helium star lifetimes as 
given by [Woosley et al. (1995 ) and compare these with 
the results obtained with the same mass-loss rates by 
Woosley et al. (1995| ) and [Wellstein fc Langer (1999| ). For 
the dashed line the final mass is obtained by numerical 
integration of the mass evolution. The final mas ses do not 
completely agree with the masses obtained by Wellstein 



Langer (1999 ), probably because these high mass-loss 



rates lead to even longer lifetimes. 

We now calculated the final masses for the revised 
mass- loss rate given by Eq. (Q), which yields 



Mf 



M 



-1.87 



2.6 X 10~*r] 



8 n -1/1.87 



(7) 



Af = -1.38 X 10"* Af 



f2.87 



(4) 



fforont from the one obtained by Nugic fe Lamoro 



and is obtained by a 'rectangular least square fit' (Langer 
1989a|]] to the data (i.e. minimising the rectangular dis- 
tances to the line, rather than the vertic al distances). The 
fit is d 
(200(f 
stars. 

For a mass-loss rate of the form 

M = -fcikf" 



In Fi g, ^we show these ma sses for the helium star lifetimes 
from Woosley et al. (1995, solid line). The lifetime of the 
helium star depends on the assumed mass-loss rate be- 
cause mass-losing helium stars become less massive, thus 
less luminous and can live longer. For example the he- 
lium star lifetimes as given by Woosley et al. (1995[) are 



substa ntially longer than the ones collected by Pols et al 



bocauDO wo o:;ccludod the hydrogen rich Wolf Rayot 



1991) for models without mass loss. We thus expect the 



(5) 



lifetimes for the helium stars with reduced mass-loss rates 
to be shorter. In Fig. Q we also plotted the final helium star 
masses assuming a lifetime which is halfway in-between 
the lifetimes given by Woosley et al. (1995 ) and Pols et al 
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Fig. 4. Final helium star masses as function of initial he- 
lium star mass with the mass-loss rates g iven by Eg. ([i| 
using a helium-star lifetime as given by Woosley ct al 
(lOOSTpolid line) and the given by Pols et al. (1991, dash- 
dotted line) and one with a helium-star lifetime halfway 
in-between these two (dashed line) . Z AMS masses for the 
helium star progenitors are indicated at the top. 



(1991, dashed line) and the one given by Pols et al. (1991 
dash-dotted line). 

The horizontal line is at 7 Mq, the typical observed 
mass of the black holes in the low-mass X-ray binaries. 
The limiting ZAMS mass for which the final helium star 
mass exceeds 7 MQwith the masses loss rate used here is 
~3O-37M0. 

We thus conclude that with revised mass-loss rates he- 
lium stars end their lives with significantly higher masses 
than previously found and may be able to form black holes 
even after case B mass transfer. 

4. Discussion 

The analysis in Sect. H neglects the influence of the wind 
of the massive star on the companion star. The compan- 
ion moves through the wind and already feels friction, 
which counteracts the widening of the orbit due to the 
stellar wind. However, even for a wind mass-loss rate of 
10~'^Mq yr^^, the density in the wind at the companion 
is almost four orders of magnitude lower than the density 
at the edge of the giant for a giant with radius of 1000 
RQand a binary separation of 1600 Rq. 

The whole argument presented in Sect. ^ is based on 
the observed mass-loss rates. However, it should be noted 
that all mass-loss rates proposed for Wolf-Rayet stars and 
used in evolutionary calculations are based on the ob- 
served rates. The valid question still remains what the 
uncertainty is in the observed mass-loss rates and in the 
inferred stellar masses and how this could influence our 
main conclusion. 



Th e mass-loss rates as determined by |Nugis fc Lamers 
(2000| ) are the most accurate, but still suffer from the gen 



& Maeder (1992) to obtain the final mass estimates from 
the luminosity. Using a different mass - luminosity rela- 
tion may change the resulting mass/mass-loss rate combi- 
nations. 

Taking the masses and mass-loss rates as plotted in 
Fig. one would not say that there is a unique mass- 
loss rate - mass relation, as is expected on theoretical 



grounds (Langer 1989b). The scatter is larger than the 
quoted uncertainty in the observations. This either points 
to underestimates of the errors in the observations, to vari- 
ability or to additional physical processes, which were not 
taken into account in the calculations by Langer (1989b ) 
and can change the mass-loss rate for a given Wolf-Rayet 
star mass. One could think of rotation, magnetic fields or 
maybe the evolutionary history. 

In the last respect it might be that stars in bina- 
ries that lose their hydrogen envelopes by mass transfer 
evolve differently from stars that lose their envelopes due 
to their own stellar winds (which possibly is enhanced by a 
companion). The question which stars actually form black 
holes and which neutron stars is considerably more com- 
plex than the question of the final mass of helium stars 
(e.g. Fryer 1999). In particular the evolution of the core 
is important. As long as the collapse of the core are not 
understood this question will remain unanswered. 

Finally, it should be noted that to form a black hole 
low-mass X-ray binary the companion must survive the 
common-envelope phase. The outcome of the common en- 
velope depends on the binding energy and density struc- 
ture of the giants envelope, which are quite different for 
giants that undergo case B and case C mass transfer. It 
could for instance be that that all binaries that undergo 
case B mass transfer to a low-mass companion will com- 
pletely merge. That would mean that we need the small 
allowed initial separation range for case C. 

5. Conclusion 



We calculated the possible initial separations for which 
case C mass transfer is likely to occur for binaries con- 
taining a massive star and a low-mass star, using different 
stellar evolution models. We find that case C mass trans- 
fer becomes impossible for primaries more massive than 
around 19 M^for th e models of Hcger et al. (2000) and 
Hurley et al. (20()0| ) and more massive than around 25 



MQfor the models by Schaller et al. (199^). For such bina- 



ries either case B mass transfer occurs, or no mass transfer 
at all. Unless the current models for massive stars underes- 
timate the radius expansion after the end of core helium 
burning the chances for forming black holes in binaries 
through case C mass transfer are therefore limited. 

We also investigated the inffuence of the assumed 
mass-loss rate on the final mass of helium stars in bina- 
ries and conclude that with a downward revised mass-loss 



rate as suggested by the observations (e.g. |Nugis fc Lamers 
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2000) helium stars end their lives with significantly higher Woosley, S. E., Langer, N., & Weaver, T. A. 1995, ApJ, 
masses than previously found and may be able to form 448, 315 
black holes even after case B mass transfer for primaries 
more massive than ^^ 30 - 40Mq. 
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